Conformally symmetric circle packings. 
A generalization of Doyle spirals 



Alexander I. BobenkoQ and Tim Hoffmann^ 

Fachbereich Mathematik, Technische Universitat Berlin, 
Strafie des 17 Juni 136, 10623 Berlin, Germany 



1 Introduction 

Circle packings (and more generally patterns) as discrete analogs of conformal mappings is a fast 
developing field of research on the border of analysis and geometry. Recent progress was initiated 
by Thurston's idea about the approximation of the Riemann mapping by circle packings. The 
corresponding convergence was proven by Rodin and Sullivan |RS]; many additional connections 



with analytic functions, such as the discrete maximum principle and Schwarz's lemma the 



discrete uniformization theorem [B£], etc., have emerged since then. 

The topic "circle packings" is also a natural one for computer experimentation and visualiza- 
tion. Computer experiments demonstrate a surprisingly close analogy of the classical theory in the 



emerging "discrete analytic function theory" |DS|. Although computer experiments give convincing 



evidence for the existence of discrete analogs of many standard holomorphic functions, the Doyle 
spirals (which are discrete analogs of the exponential function, see section ^ are the only circle 
packings described explicitly. 

Circle packings are usually described analytically in the Euclidean setting, i.e. through their 
radii function. On the other hand, circles and the tangencies are preserved by the fractional-linear 
transformations of the Riemann sphere (Mobius transformations). It is natural to study circle 
packings in this setting, i.e. modulo the group of the Mobius transformations. He and Schramm 
p3S(| developed a conformal description of hexagonal circle packings, which helped them to show 
that Thurston's convergence of hexagonal circle packings to the Riemann mapping is actually . 
They describe circle packings in terms of the cross-ratios 

{a-h){c-d) 

q{a,b,c,d) := — — 

(o — c){a — a) 

of their touching points. 

In Schramm introduced circle patterns with the combinatorics of the square grid (SG pat- 
terns). In many aspects the SG theory is analogous to the theory of the hexagonal circle packings. 
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Figure 1: A circle flower. 



However, the SG theory is analytically simpler. The corresponding discrete equations describing 
the SG patterns, in the Euclidean as well as in the conformal setting, turn out to be integrable 
I BP]. Methods of the theory of integrable equations allowed us in ^ | |AB| ] to find Schramm's circle 



patterns which are analogs of the holomorphic functions 2;",logz. 

One big question is which results on the Schramm's circle patterns carry over to the hexagonal 
setting, in particular whether some discrete standard functions can be described explicitly. This is 
closely related to the question of integrability of the basic discrete equations for hexagonal circle 
packings (He-Schramm equation, see section ^). In the present paper the first simple step in this 
direction is made. We study (surprisingly non-trivial) conformal geometry of hexagonal circle 
packings. In terms of this approach, a special class of conformally symmetric circle packings, which 
are generalizations of the Doyle spirals, is introduced and all such packings are described explicitly. 

Since this paper deals with families of circle packings it seems natural to show not only arbi- 
trarily choosen members in the figures, but to provide a possibility to present them all. Therefore 
there is an interactive version of this paper available [online]. It has some of the figures replaced 
by applets, that allow to explore the families directly. See section ^ for more information on this 
version. 

Acknowledgments. The authors thank U. Hertrich-Jeromin, U. Pinkall, Yu. Suris and E. Tjaden 
for helpful discussions. 



'Discrete z'^ and log z have been conjectured by Schramm and Kenyon earher [ 3www| 
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2 Geometry of circle flowers and conformaUy symmetric circle 
packings 

This paper concerns patterns of circles in the plane called hexagonal circle packings. Their basic 
unit is the flower, consisting of a center circle tangent to and surrounded by petals. A hexagonal 
flower is illustrated in Fig|l|; the six petals form a closed chain which wraps once in the positive 
direction about the center. Whereas the neighboring petals touch, the circles of not-neighboring 
petals of a flower may intersect. We call a flower immersed if none of its circles degenerates to a 
point. A hexagonal circle packing is a collection of oriented circles where each of its internal circles 
is the center of a hexagonal flower. Orientations of the circles should agree: at the touching points 
the orientations of the touching circles must be opposite. A hexagonal circle packing can be labeled 
by the triangular (hexagonal) lattice 

HL = n + me'''/^ G C, n,m£Z. 

or by one of its subset. A circle packing is called immersed if all its flowers are immersed. Immersions 
of the whole HL are called entire. Fractional-linear transformations of the complex plane (Mobius 
transformations) preserve circles, their orientation and their tangencies. In this paper we study 
circle packings modulo the group of Mobius transformations. 

The center circle of a flower contains 6 points zi, . . . , zq (see Fig.|l]) where it touches the petals. 
We call them center touching points of a circle flower. 

Proposition 2.1 Let z\,... , zq be cyclicly ordered^ points on a circle C. Then the following three 
statements are equivalent: 

(i) There exists a flower with the center C and center touching points zi, . . . ,zq, 

(a) The multi-ratio m of zi, . . . ,zq is equal to —1, i.e. 

m{zi,Z2, zs, Z4, Z5, zq) := = -1, (1) 

{Z2 - Z3){Z4 - 2:5) (^6 - Zl) 

(Hi) There exists an involutive Mobius transformation M (Mobius involution) such that 

M{zk) = Zk+3 (k mod 6). 



Proof. Mapping the point zq to infinity by a Mobius transformation one obtains two parallel lines 
and five touching circles as in Fig.|2[ An elementary computation yields 

Zk+i - Zk = 2^rk+irk, A; = 1, . . . , 4, (2) 

where are the radii of the corresponding circles. Together with ri = r^ and (z^—zq) / {zq—zi) = —1 
this implies (||). 

^Positive orientation of the ordering is assumed. 
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Zi Z2 Z3 Z5 

Figure 2: A flower with one central touching point at infinity. 



On the other hand, given arbitrary ri > and ordered zi, . . . ,zq satisfying (|l|), after normahzing 
zq = 00 formula (|2|) provides us with the radii of the touching circles as in Figj2|. This proves the 
equivalence of (i) and (ii). 

To show the equivalence of (ii) and (iii), define the Mobius transformation M through M{zi) = 
Zi, M{z2) = 2:5, M{z^) = zq. Consider = M{zi). The invariance of the cross-ratios q{zi, Z2, z^, Z4) 
= q{z4, z^, zq, z^) implies the equivalence of (||) and = zi. The same proof holds for M(z5) = Z2 
and M{zq) = z\. 

To each center touching point of a fiower, one can associate a circle passing through 4 
touching points Zfc-i, z^+i, u)fc, of the fiower containing z^ (see FigJ^. Here is the touching 
point of petals^ "Pfc+i and V^. Indeed, mapping the point by a Mobius transformation to cxd, it 
is easy to see that the points Zfc_i, z^+i, uifc, are mapped to vertices of a rectangle, thus lie on 
a circle. We call these circles s-circles of a fiower. 

Theorem 2.2 There exist a one-parameter family of flowers with the same center touching points. 
Moreover, there exists a unique flower F in this family, which satisfies the following equivalent 
conditions: 

(i) F is invariant with respect to a Mobius involution M with a fixed point P, 
(ii) All s-circles of F intersect in one point P. 

We call the fiower F of the theorem conformally symmetric. 
One can view the whole family of fiowers at [appl] . 

Proof. Keeping the points z\,. . . , Z5 in Fig.^ fixed and varying r\ one obtains a one param- 
eter family of flowers with the same touching central points. Let us now construct the fiower 



F . The Mobius involution of Proposition 2.1 preserves the central circle C. Consider the circles 



The petals are labeled by the corresponding touching points Zh- 
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Figure 3: A conformally symmetric flower. 



Cfc, i = 1, 2, 3, orthogonal to C and passing through the pairs of points {z^, -Zfc+s}. All these three 
circles intersect in 2 points P and P', which are the fixed points of M lying inside and outside C, 
respectively. By a Mobius transformation, map the point P to infinity. The Mobius involution M 
becomes M{z) = —z and the circles Ci,C2,C3 become straight lines intersecting in the center of 
C . To construct the flower connect the z^-points with even (respectively, with odd) labels by 
straight lines and consider their intersection points (see Fig.^). The circles Cfc passing through 
the triples Wh^Wh-\^z\^ touch at the points w^. Let us prove this fact for C\ and C^- Indeed, the 
triangles ^{wx^w^^zx) and A (23, 2:5, zi) are similar, therefore the tangent lines to the circle C\ at 
w\ and to the circle C at Z3 are parallel. The tangent lines to C2 at w\ and to C at zg are also 
parallel. Since the points z^ and zg are opposite on C, the circles C\ and C2 touch at w\. The 
circles are the petals of the desired flower which is obviously M-symmetric. The s-circles of 
this flower are the straight lines (^^,2^+2)- The latter obviously intersect at infinity, thus all the 
s-circles of F intersect in the fixed point P of M. 

The proof of {ii) =^ (i) is similar. After mapping the point P to infinity the s-circles become 
straight lines and the flower is as in Fig.^ Since the circles in this figure touch, their tangent lines 
at the points Zfc,Zfc+3 and tf^fc+i are parallel. This implies that z^ and 2^+3 are opposite points on 
C, and the flower is symmetric with respect to the 7r-rotation of C . 

Definition 2.3 A hexagonal circle 'packing is called conformally symmetric or an s-circle packing 
if it consists of conformally symmetric flowers, i.e. the s-circles of each of its flowers intersect in 
one point. 
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Figure 4: A normalized confer mally symmetric flower. 



3 Analytic description of conformally symmetric circle packings 

In this section we describe all conformally symmetric circle packings using the conformal description 
of circle packings proposed by He and Schramm |HS|. 

To each central touching point of a flower one associates the cross-ratio^ 

, ^ {Zk - Zk-l){Wk-l - Wk) 

[Zk-l - Wk-lj{Wk - Zk) 

Mapping to oo, one observes that tree other points in @ are mapped to vertices of a rectangle, 
which implies that Sk is purely imaginary. Moreover, the cross-ratios of an immersed oriented flower 
are positive imaginary, —isk > 0. Also note that 

Sk = -q{zk+i,Zk-i,Wk-i,Zk) = q(zk,Wk,Zk+i,Zk^i), (4) 

and that = q{zk+i, Zk-i,Wk-i,Wk) is the cross-ratio of the four touching points lying on the 
s-circle Sk- 



Lemma 3.1 The cross-ratios Sk of a flower satisfy the He-Schramm equation [HS] 

Sk + Sk+2 + Sk+4 + SkSk+lSk+2 = (5) 



^Note that our normalization of Sk differs from the one in [HS]. 
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for all k mod 6. 



Proof. Let be the Mobius transformation that takes z^, z^^i, Wk^i to the points oo, 0, 1, respec- 
tively. By the definition of Sfc we have 

Sk = qizk,Zk^i,Wk-i,Wk) = q{oo,0,l,mk{wk)), 
-Sk = q{zk+i,Zk-i,Wk-^i,Zk) = g(mfc(zfc+i),0, l,oo), 

thus 

mkiwk) = 1 - Sk, mk{zk+i) = -Sk- 
For Mk := nik+imi^ this yields Mk{—Sk) = oo, Mfc(oo) = 0, Mk{l — Sk) = 1 and, finally, 

where the usual matrix notation for the Mobius transformation is used. The equality of the corre- 
sponding Mobius transformations implies M3M2M1 = ±M^^ M^^ M^^ , which is 

/ S2 1 + S1S2 ^ — ± - 56 - S4S5S6 1 + S4S5 \ 

\ 1 + S2S3 S1 + S3 + S1S2S3 J \ 1 + S5S6 -S5 J ■ 

Since the set of immersed flowers is connected and s's do not vanish the sign in this equation is the 
same for all flowers. Taking all the circles with the same radius one checks that the correct sign is 
plus, which implies the claim. 

Given a hexagonal circle pattern it is convenient to associate its touching points as well as the 
cross-ratios to the edges of the honeycomb lattice. Equation (|5|) is a partial difference equation 
on the honeycomb lattice. The cross-ratios on the edges of each hexagon satisfy (|5|). Moreover, 
it is easy to check that the He-Schramm equation is sufficient to guarantee the existence of the 
corresponding circle packing. 

Proposition 3.2 Given a positive-imaginary function s : E ^ iM+ on the edges E of the honey- 
comb lattice satisfying ^ on each honeycomb, there exists unique (up to Mobius transformation) 
immersed hexagonal circle packing with the cross-ratios given by the corresponding values of s. 

Theorem 3.3 A circle flower is conformally symmetric if and only if its opposite cross-ratios s^ 
are equal 

Sk = Sfc+3 {k mod 6). (6) 

Proof. The property (|6|) for conformally symmetric flowers follows from (i) of Theorem A 
simple computation with the flowers in FigJJ shows that the map (si, S2) of immersed conformally 
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symmetric flowers to (iM+)^ 9 {81,82) is surjective. Since a flower is determined through the s's, 
the converse statement follows. 

The general solution of (P, |6|) on the whole HL depends on three arbitrary constants and can 
be given explicitly. There is a JAVA applet that lets you explore this three parameter family of 
circle packings interactively at [appl]. 




Figure 5: Cross-ratios of conformally symmetric circle patti 



Theorem 3.4 The general solution of ^) is given by 

Qn = itan(An + a), 

bn = itan{An + (3), (7) 
Cn = itan(An + 7), 

where A = —a — /? — 7 and the cross-ratios Sk on the edges of the hexagonal lattice are labeled by 
an,bn,Cn as shown in Fig.^. 

Proof. We start with a simple proof of the consistency of the following continuation of a solution 
of P). Given s satisfying on a honeycomb H, i.e. a, b, a in Fig.|6| satisfying 

a + b + c + abc = 0, (8) 

and a value of s on one of the edges attached to the honeycomb ( for example, di in Fig.E), it can 
be uniquely extended to the full six honeycombs Hi, . . . ,Hq neighboring H. Indeed, (P, t) yield 

b + di+d2 + bdid2 = 0, 
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Figure 6: Continuation of conformally symmetric s about a honeycomb. 



thus d2 = Mi{di) is a Mobius transformation of di. Passing once around the honeycomb H in this 
way one can check that (|8|) imphes the monodromy Mobius transformation M = Mq . . . Mi is the 
identity, thus this continuation implies no constraints on di. 

Proceeding this way, one reconstructs s on the whole lattice HL from its values on three adjacent 
edges {a,b,di above). Then ^) imply 

Qn + b_n + Ci + Qnb-nCl = 0, Qn+l + b-n + Cq + a„+i6_„Co = 

and similar relations for other an,bn,Cn. These identities become just the addition theorem for the 
tangent function, implying the formulas in (^), which can be checked directly. 



4 The Doyle spirals 

Denote by R the radius of the center circle of a flower and hy Rk, k = 1, . . . ,6, the radii of its 
petals. The Doyle spirals are characterized through the constraint (see | BDS| , CR] for a complete 
analysis of Doyle spirals) 

RkRk+3 = R^-. RkRk+2Rk+4 = (9) 

on the radii of the circles (see Fig.0 where the central radius is normalized to be i? = 1). The Doyle 
spirals have two degrees of freedom (for example the ratios Ri/R and R2/R, which are the same 
for the whole spiral) up to similarities. Again, you can experiment with the two radii in a JAVA 
applet [appl]. 

Proposition 4.1 The Doyle spirals are conformally symmetric. 

Proof. The configurations of four touching circles with the radii R, Rk-i, Rk, Rk+i stnd with the 
radii Rk+3, Rk+i, R, Rk+2 differ by scaling. This implies Sk = Sk+3 (use both and the second 
representation of Sk in (^)) and the claim follows by Theorem 3.3 . 
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Figure 7: Radii of a Doyle spiral with the normalized central radius R = 1. 



Theorem 4.2 The Doyle spirals and their Mobius transformations can be characterized by the 
following two equivalent properties: 

(i) The circle packing is conformally symmetric, and the corresponding solution of is "con- 
stant". It is of the form ^ with a,P,j S (0, a + /? + 7 = (mod tt) or, equivalently, 

a„ = ao, bn = bo, c„ = cq, oq, 60, cq G ao + bo + cq + aoboco = 0. 

(a) The whole circle packing is invariant with respect to the Mobius involution of each of its 
flowers. 

Proof. All the flowers of a Doyle spiral differ by scaling, which implies (i). Consider the Doyle 
spiral as in FigJ^. Computing the cross-ratios through the radii, one shows that the map 

{{A, B) G M^} ^ {(a, 6, c) G : a + b + c + abc = {)}. 

is surjective, thus (i) characterizes the Doyle spirals and their Mobius transforms. The proof of the 
equivalence (i) 4^ [ii) is elementary and is left to the reader. 

It is an open problem whether the Doyle spirals are the only entire circle packings. Formulas 
imply that it is possible to have all cross-ratios being positive imaginary (necessary condition 
for entireness) only when A = 0. 

Corollary 4.3 Doyle spirals are the only entire conformally symmetric circle patterns. 



10 



5 Airy functions as continuous limit 

Because of the property (^), Doyle spirals are interpreted as a discrete exponential function. 

In the conformal setting this interpretation can also be easily observed. Indeed, let be a 



family of circle packings approximating a holomorphic mapping in the limit e ^ 0. In [HS] He and 
Schramm investigated the behavior of the cross-ratios Sk in this limit: 



Sk 



where h^. is called the discrete Schwarzian derivative (Schwarzian) of P*" at the corresponding edge 
of the hexagonal lattice. The discrete Schwarzians converge to the Schwarzian derivative 

of the corresponding holomorphic mapping. More precisely, there exist continuous limits 

a = lim /ii , b = lim h%, c= lim h\ 

for the smooth function^ a, b, c. Because of (^) these functions satisfy 

a + 6 + c = (11) 

at each point. The Schwarzian equals 

S{f)=i{a + qH + qc), q = e^^'l\ 



and, using (|ll|) , this also yields 

6a = Re(5(/)), Qb = Re{qS{f)), Qc = Re{q^S{f)). (12) 



We see that, due to Theorem f4.2| , Doyle spirals correspond to holomorphic functions with constant 
Schwarzian derivative S{f) = const. The general solution of the last equation is the exponential 
function and its Mobius transformations. 

It is natural to ask which holomorphic functions correspond to general conformally symmetric 
circle packings. In Figjs] one observes that each of the cross-ratios a„,6„,c„ is constant along one 
lattice direction. For the functions a, 6, c above, this implies 

a = a(Re(z)), b = &(Re(gz)), c = c{Yle{q^z)), (13) 

where z is the complex coordinate. Comparing (|l^ ) and (^) we see that the Schwarzian is a linear 
function of z: 

S{f) = Az + B, AeR, B eC. (14) 



^Note that lim(,_,o hi — lirriE^o ^fc+a 
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Equation (|T^ can be easily solved by standard methods. The general solution of S{f) = u{z) with 
holomorphic u{z) is given by f{z) := V'i/V'2, where Tpi{z) and ip2{z) are two independent solutions 
of the linear differential equation ip" = u{z)'iIj. 

By a shift and scaling of the variable z, equation (|T4[) with A ^ can be brought to the form 



S{f) = z. (15) 
Solutions of the corresponding linear equation 

ip" = zip (16) 



are the Airy functions Ai(2;) and Bi(2;). On the real line the first one is given by |S0] 

Ai(x) = — / cos{xt H )dt, 

Jo 3 

and the second one is related to it by 

Bi(2:) = iq^A\{q^z) - iqA\{qz). 
In the corresponding Mobius class of solutions of ([l5|) it is natural to choose 



which is the most symmetric one, f{qz) = qf{z). The corresponding circle packing, symmetric with 
respect to the rotation z qz, is shown in Fig.^. 

6 About the interactive version of this paper 

Since the families of circle packings discussed in this paper have a finite (and even small) number 
of parameters it seemed to be natural to look for a way to visualize the whole families. 

Except the Doyle spirals the families are only defined modulo an arbitrary Mobius transforma- 
tion. Therefore it should be possible to view the packings with a Mobius transformation applied, 
so one could look at them through "Mobius glasses" . 

The outcome of this is an interactive version of this paper. This version includes some Java 
applets that let you experiment with the circle packings directly. In particular there are applets 
to illustrate the families of circle packing flowers, the whole class of conformally symmetric circle 
packings, and the special case of Doyle spirals. 

This interactive version renders a dvi file inside a Java applet. It needs a web browser that 
includes a Java vm. Since flipping pages might get slow on old machines and Java vm's that have 
no just in time compiler, we also provide a page that shows the applets that are missing in this 
version only. You will find the paper and the applets at 



http://www-sfb288.math.tu-berlin.de/Publications/online/ 
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Figure 8: A conformally symmetric circle packing (with a = /3 = 7 in (|^)) and its smooth counter- 
part. The vertices of the hexagons are the images of the points of a standard hexagonal grid under 
the map / from ([TtI). 
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